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Abstract 

Embedding diagrams prove to be quite useful when learning general relativity as they offer a 
way of visualizing spacetime curvature through warped two dimensional (2D) surfaces. In this 
manuscript we present a different 2D construct that also serves as a useful conceptual tool for 
gaining insight into gravitation, in particular, orbital dynamics - namely the cylindrically sym¬ 
metric surfaces that generate Newtonian and general relativistic orbits with small eccentricities. 
Although we first show that no such surface exists that can exactly reproduce the arbitrary bound 
orbits of Newtonian gravitation or of general relativity (or, more generally, of any spherically sym¬ 
metric potential), surfaces do exist that closely approximate the resulting orbital motion for small 
eccentricities; exactly the regime that describes the motion of the solar system planets. These 
surfaces help to illustrate the similarities, as well as the differences, between the two theories of 
gravitation (i.e. stationary elliptical orbits in Newtonian gravitation and processing elliptical-like 
orbits in general relativity) and offer, in this age of 3D printing, an opportunity for students and 
instructors to experimentally explore the predictions made by each. 
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I. INTRODUCTION 


One hundred years ago, Albert Einstein published his general theory of relativity (gen¬ 
eral relativity or GR), which describes gravity as the warping of space and time due to the 
presence of matter and energy, and forever changed our conceptual understanding of one 
of nature’s fundamental forces. Prior to GR, Isaac Newton’s theory of gravitation, which 
describes gravity as the force of one massive object on another, reigned supreme and unchal¬ 
lenged for nearly a quarter of a millennium. With the advent of GR, Newtonian gravitation 
was reduced to that of an approximation to a more fundamental underlying theory of na¬ 
ture and must be abandoned altogether when probing phenomena near highly relativistic 
objects such as black holes and neutron stars. Nevertheless, Newtonian gravitation is highly 
successful in accounting for the approximate motion of the planets in our solar system and 
allows for a theoretical derivation of Kepler’s three empirical laws. 

According to Newtonian gravitation (and Kepler’s first law of planetary motion) the 
planets move in closed elliptical orbits with the Sun residing at one focus. [1] For these orbits, 
the translating body hnds itself at precisely the same radial distance after the polar angle has 
advanced by 27r and will therefore retrace its preceding elliptical motion indefinitely. This 
Newtonian prediction of closed elliptical orbits holds true only in the absence of other massive 
objects as additional orbiting bodies generate perturbations to the motion. Additionally, 
when the theoretical framework of GR is employed, the closed or stationary elliptical orbits 
of Newtonian gravitation are replaced with elliptical-like orbits whose perihelia and aphelia 
precess, even in the absence of additional massive objects. [2] Before Einstein’s GR, the small 
but observable precession of Mercury’s perihelion was partially unaccounted for and offered 
the hrst real test of general relativity. In fact, when the precession of Earth’s equinoxes 
and the gravitational effects of the other planets tugging on Mercury are accounted for, 
GR claims exactly the additional measured precession. The predicted value of Mercury’s 
perihelion precession due to general relativistic effects is incredibly tiny, equating to an 
angular shift per orbit of less than one part in ten million. 

According to general relativity, a massive object causes the spacetime in its vicinity to 
warp or curve simply due to its presence. Hence, the planets move in their respective orbits 
not due to the gravitational force of the Sun acting on them, but rather due to the fact 
that the planets are moving in the warped four-dimensional spacetime about the Sun. In 
a typical introductory undergraduate GR course, the spacetime external to a non-rotating, 
spherically symmetric massive object is often extensively studied as it offers a simple, exact 
solution to the held equations of general relativity. To gain conceptual insight into this 
non-Euclidean spacetime, an embedding diagram is often constructed. By examining the 
spacetime geometry external to the massive object at one moment in time, with the spherical 
polar angle set to 7r/2, a two dimensional (2D) equatorial “slice” of the three dimensional 
space is extracted. A curved two dimensional, cylindrically symmetric surface embedded in 
a hat three dimensional space is then constructed with the same spatial geometry as this 
two dimensional spacetime slice; this is the embedding diagram. These embedding diagrams 
prove incredibly useful in visualizing spacetime curvature and serve as a conceptual tool. 

When free particle orbits, or geodesics, about non-rotating, spherically symmetric massive 
object are then later studied, a conceptual analogy of a marble rolling on a warped two- 
dimensional surface is often drawn to represent the orbital motion. Although this analogy 
may be useful for the beginning student of GR in visualizing particle orbits about massive 
objects in curved spacetime, it isn’t precise and fundamentally dihers from the analogy 
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of embedding diagrams discussed above. Although the function that describes the shape 
of the cylindrically symmetric surface of the embedding diagram is of a relatively simple 
form, there exists no 2D cylindrically symmetric surface residing in a uniform gravitational 
field that can generate the exact Newtonian orbits of planetary motion, save the special 
case of circular orbits. [3] This null result holds when considering a fully general relativistic 
treatment of orbital motion about a non-rotating spherically symmetric massive object. [4] 
For an object rolling or sliding on the interior of an inverted cone with a slope of h has 
been shown that the stationary elliptical orbits of planetary motion can occur when the orbits 
are nearly circular. [5] Later, elliptical-like orbits with small eccentricities were considered 
on 2D cylindrically symmetric surfaces where a perturbative solution was constructed. [6] 
There it was found that when the surface’s shape takes the form of a power law, precessing 
elliptical-like orbits with small eccentricities can only occur for certain powers and stationary 
elliptical orbits with small eccentricities can only occur for certain powers and for certain 
radii on the surface. 

Although we show that no such surface exists that can exactly reproduce the arbitrary 
bound orbits of a spherically-symmetric potential, cylindrically symmetric surfaces do ex¬ 
ist that yield the stationary and precessing elliptical-like orbits with small eccentricities of 
Newtonian gravitation and general relativity, respectively. First, we arrive at a general ex¬ 
pression describing the slope of a 2D cylindrically symmetric surface that yields the elliptical 
orbits of Newtonian gravitation for small eccentricities. The slope of this surface reduces to 
the special case of White’s when the arbitrary constant of integration is set equal to zero. [5] 
We then extend our study to general relativity and arrive at an expression for the slope of 
a 2D surface that generates the precessing elliptical orbits of GR for small eccentricities. 
For the sake of brevity we symbolically represent stationary elliptical-like orbits of Newto¬ 
nian gravitation with small eccentricities as Ne orbits and precessing elliptical-like orbits of 
general relativity with small eccentricities as GRe orbits; these phrases appear frequently 
throughout the paper and warrant this use of shorthand. 

This paper is outlined as follows. In Sec. II, we derive the equations of motion for an 
object orbiting in a generic spherically symmetric potential. We examine the case of Newto¬ 
nian gravitation and present the solution for the respective particle orbits while introducing 
the notation contained in the sections that follow. In Sec. Ill, we consider the rolling or 
sliding motion of an object constrained to reside on a 2D cylindrically symmetric surface 
and obtain the corresponding orbital equation of motion. We show that this equation is 
fundamentally different than the equation of motion for an object subjected to a generic 3D 
spherically symmetric potential and that there exists no 2D cylindrically symmetric surface 
that can make these equations coincide for any spherically symmetric potential, except for 
the special case of circular orbits. For elliptical-like orbits with small eccentricities on the 
cylindrically symmetric surface, we solve the orbital equation of motion perturbatively. A 
differential equation, evaluated at the characteristic radius of the orbital motion, emerges 
and relates the precession parameter of the orbit to the slope of the respective surface. By 
demanding that the elliptical-like orbits on the surface have a constant precession parameter 
and by relaxing the condition that the aforementioned differential equation be evaluated at 
the characteristic radius of the motion, we solve the corresponding differential equation and 
find the slope of the surface that yields these elliptical-like orbits with small eccentricities 
for all radii on the surface. We then examine the solution for the slope that generates the 
N5 orbits. In Sec. IV, we extend our study to general relativity. By demanding that the 
elliptical-like orbits on the 2D surface now mimic those of GR, we present the slope of the 
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surface that yields these GRe orbits for all radii on the surface. We then compare these 
two families of surfaces with the intent of offering insight into the aforementioned theories 
of gravitation. 


II. THE CLASSICAL EQUATIONS OF MOTION FOR AN OBJECT ORBITING 
IN A SPHERICALLY SYMMETRIC POTENTIAL AND THE SOLUTION FOR A 
NEWTONIAN POTENTIAL 


We begin this manuscript with a review of central-force motion and, in particular, of 
Newtonian gravitation. This treatment is not only instructive but also serves to introduce 
the notation contained within the sections that follow. 

The Lagrangian describing an orbiting body of mass m in a generic spherically symmetric 
potential is of the form 


L 


1 

-m 


sim 


U(r), 


( 2 . 1 ) 


where r,9,(j) are spherical polar coordinates and a dot indicates a time-derivative. As the 
angular momentum of an orbiting body in a spherically symmetric potential is conserved, 
the motion of the object, withont loss of generality, can be chosen to reside in the equatorial 
plane. This is done by setting 6 = 7r/2 in Eq. (2.1) where the Lagrangian rednces to one 
dependent only on the two plane polar coordinates r and 0. It is noted that conservation of 
angnlar momentnm is a general result of central-force motion and, for the case of planetary 
motion in Newtonian gravitation, eqnates to Kepler’s second law. [7] 

The Lagrangian of Eq. (2.1) generates two eqnations of motion, which take the form 

^2 1 dU{r) , , 

^ ^ ^ ( 2 . 2 ) 


r -5- + 


m dr 


= 0 


4 ' = L (2.3) 

where £ is a constant of the motion and eqnates to the conserved angular momentum per 
unit mass. Here we are interested in arriving at an expression for the radial distance of the 
orbiting body in terms of the azimuthal angle, r{(j)). Using the chain rnle and Eq. (2.3), a 
differential operator of the form 


d id 
dt dcj) 

can be constructed and used to transform Eq. (2.2) into the form 

d‘^r 2 i dr'\^ r^ dU{r) 

dcj)^ r ^ ^ dr 


(2.4) 


(2.5) 


Equation (2.5) is the eqnation of motion describing the radial distance of an orbiting body 
as a fnnction of the azimnthal angle in a generic spherically symmetric potential. This 
eqnation will later be compared to the eqnation of motion for an object rolling or sliding on 
a cylindrically symmetric 2D snrface. 

For an object of mass m residing in a Newtonian gravitational potential, the general 
expression of Eq. (2.5) takes the form 


^ _ 2 2 _ Q 

dcj)^ r P 


( 2 . 6 ) 
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where G is Newton’s constant and M is the mass of the central object. Equation (2.6) is 
usually presented in a slightly simpler form in most classical dynamics texts by adopting 
a change of variables of the form u = 1/r. This substitution is unnecessary in the present 
treatment and will thus be avoided. It is noted that Eqs. (2.2) and (2.3) are technically 
the equations of motion for the reduced mass /i of the equivalent one-body problem, but 
approximately describe the motion of the particle of mass m in the limit that M 3> m.[7] 
In this limit, we can neglect the motion of the central mass M treating it as hxed, which is 
exactly the regime this manuscript is interested in understanding. 

Equation (2.6) can be solved exactly and yields the well-known conic sections of the form 


r(0) 


^0 

(1 -|- e COS0) ’ 


(2.7) 


where tq is a constant, known historically as the semilatus rectum, and e is the eccentricity 
of the orbit. Notice that for the eccentricity, £ = 0 corresponds to the special case of circular 
motion whereas 0 < e < 1 describes elliptical orbits. Also notice that Eq. (2.7) is a solution 
to Eq. (2.6) when = GMtq. This relation specihes the necessary angular momentum at 
a given radius, characterized by tq, for circular or elliptical motion to occur. 

For small eccentricities, the exact solution of Eq. (2.7) can be approximated as 


r{(l>)ap = ro{l-e cos 0) (2.8) 

to hrst-order in the eccentricity. Notice that when Eq. (2.8) is taken as the solution for the 
radius, tq corresponds to the average radius of an elliptical orbit with small eccentricity to 
hrst-order in the eccentricity. 

The eccentricities of the planets of our solar system are quite small with an average 
eccentricity of Save = 0.060, hence Eq. (2.8) yields an excellent approximation for the radii 
of the orbits. As a specihc example. Mercury has the largest of eccentricities with a value 
of Em = 0.2056. Comparing the exact and approximate radii given by Eqs. (2.7) and (2.8) 
respectively, we hud a maximum percent error of 4.2%.[8] See Table I near the end of this 
article for the eccentricities of the solar system planets. 


III. THE CLASSICAL EQUATIONS OF MOTION FOR AN OBJECT ORBITING 
ON A CYLINDRICALLY SYMMETRIC SURFACE 


The Lagrangian describing a body of mass m residing on a cylindrically symmetric surface 
in a uniform gravitational held is of the form 


L = -,n 


•2 I 2 12 I -2 

r + r 0 + z 


+ - mgz, 


(3.1) 


where we have included both the translational and rotational kinetic energy of the orbiting 
body and the gravitational potential energy of the object at a height 2 :. In an earlier version 
of this manuscript, the rolling and sliding motions of an orbiting body were collectively 
analyzed and presented. In the analysis the rolling object was assumed to roll without 
slipping and the approximation that the angular velocity of a rolling object is proportional 
to the translational velocity of the object’s center of mass was adopted. There it was found 
that the general solution for the slope of the surface, which will be presented later in this 
section, does not depend on the intrinsic spin of the orbiting body. In fact, none of the 
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results presented in this manuscript depend on the intrinsic spin of the orbiting body. With 
the benefit of hindsight and for the sake of simplicity, we forego this more detailed analysis 
as the results are left unaltered and set the rotational kinetic energy term of Eq. (3.1) to 
zero. 

As we are interested in studying the motion of an orbiting body on a cylindrically sym¬ 
metric surface, we note that there exists an equation of constraint connecting two of the 
three cylindrical coordinates generically of the form = z{r). Using this fact and setting 
the rotational kinetic energy term to zero, Eq. (3.1) can be massaged into the form 

L = (1-I-U(r)^)r^-I- —mgz{r), (3.2) 

where z'{r) = dz{r)/dr and we used the fact that i = rz' via the chain rule. This Lagrangian 
generates two equations of motion, which have been previously derived elsewhere and are of 
the form 

(1 + U^)r -I- z'z''r^ — ^ + gz' = Q (3.3) 

(3.4) 

where, again, ^ is the conserved angular momentum per unit mass. [3, 4] 

Upon comparison of Eqs. (2.2) and (3.3), we note that there exists no spherically sym¬ 
metric potential, U(r), whose presence exists in the form of a radial derivative in Eq. (2.2), 
that can generate the second term of Eq. (3.3). This is easily seen as the last term of 
Eq. (2.2) will only generate functions of r and not r. Likewise, there exists no cylindrically 
symmetric surface, described by z{r), that can generate a term capable of canceling the 
term in Eq. (3.3) as derivatives of z{r) are likewise strictly functions of r. The term 
of Eq. (3.3) arises from motion in the ; 2 -direction on the cylindrically symmetric surface 
where r 7 ^ 0. Hence, Eqs. (2.2) and (3.3) cannot be made to match for any spherically 
symmetric potential, U(r), or for any 2D cylindrically symmetric surface, z(r), except for 
the trivial case of U'{r) = z'{r) = 0. This disparity implies that there exists no cylindrically 
symmetric surface residing in a uniform gravitational held that is capable of reproducing the 
exact motion of a body orbiting in any spherically symmetric potential U = U{r), except 
for the special case of circular motion. This null result was hrst pointed out for an orbiting 
object subjected to the Newtonian gravitational potential of a central mass and later for 
the case of a body orbiting about a non-rotating, spherically symmetric massive object in 
general relativity. [3, 4] For the special case of circular motion, a 2D cylindrically symmetric 
surface can yield the orbits of Newtonian gravitation, taken here to mean that the surface 
yields orbits that obey Kepler’s three laws, when the shape of the surface takes the form of 
a — 1/r Newtonian gravitation potential well. This result only holds true for circular motion 
as elliptical-like orbits on this surface precess for all radii and hence violate Kepler’s hrst 
law. [ 6 ] See Fig. 2(a) for a surface of revolution plot of z{r) vs r of a Newtonian gravitational 
potential well. 

There exists no cylindrically symmetric surface residing in a uniform gravitational held 
capable of generating the exact Newtonian orbits of planetary motion, except for the case 
of circular orbits. This can be witnessed explicitly by inserting the solution for Newtonian 
orbits, given by Eq. (2.7), into Eq. (3.3). Upon this substitution, Eq. (3.3) takes the form 

3 

(1 — %^z') -\-e{2 — z'"^) cos + [(1 ~ 2z'‘^) cos^ 0 — ^0 sin^ <p z'z”] — cos^ p = 0, (3.5) 

£2 L J 
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where we employed the chain rule and Eq. (3.4) in calculating time derivatives and grouped 
the resultant expression by powers of the eccentricity. In order for Eq. (2.7) to be an 
exact solution to Eq. (3.3) for arbitrary eccentricity, each of the four terms of Eq. (3.5) 
must vanish uniquely. Although the first two terms can in fact be set equal to zero, which 
consequently would yield the required angular momentum per unit mass for elliptical or 
circular motion to occur and a slope for the surface, notice that the and terms will not 
vanish for any cylindrically symmetric surface. For circular orbits, where the eccentricity is 
set to zero, Eq. (2.7) represents an exact solution to Eq. (3.3) when the angular momentum 
per unit mass obeys the expression = gr^z' for a given characteristic radius and slope. 


A. The 2D surfaces that generate elliptical-like orbits with small eccentricities 


Using Eq. (2.4), Eq. (3.3) can be transformed into an orbital equation of motion of the 
form 

Equation (3.6) is a non-linear equation of motion that can be solved perturbatively. For 
elliptical-like orbits with small eccentricities, we choose an approximate solution for the 
radius of the form 

r(0) = fo(l - £ cos(z/0)), (3.7) 

where fg and u are parameters that are to be determined by Eq. (3.6) and e is the eccentricity 
of the orbit, which will be treated as small and used as our expansion parameter.[9] As 
the remainder of this paper deals solely with elliptical-like orbits of small eccentricity, all 
subsequent orbits should be understood as having small eccentricities, even if not explicitly 
stated. 

The precession parameter, z/, is understood as 



(3.8) 


where A 0 corresponds to the angular separation between two apocenters of the orbital 
motion. [2, 6 ] Notice that when z/ = 1, Eq. (3.7) has the same functional form as that 
of Eq. (2.8) and describes a stationary elliptical-like orbit on the 2D surface where the 
angular separation between the two apocenters is 27r. When z/ 7 ^ 1, Eq. (3.7) describes a 
processing elliptical-like orbit on the 2D surface where A0 7 ^ 2tt. Interestingly, orbits about 
non-rotating, spherically-symmetric massive objects in general relativity are found to have 
precession parameters that take on values less than one. This corresponds to an angular 
separation of A 0 > 27r, where the apocenter of the elliptical-like orbit marches forward in 
the azimuthal direction in the orbital plane. 

Inserting Eq. (3.7) into Eq. (3.6) and keeping terms up to hrst-order in e, we hnd that 
Eq. (3.7) represents an approximate solution to the orbital equation of motion when 

(3-9) 

z'^{l + z'^y = ^z'^ + rozl (3.10) 

where Zq and Zq are radial derivatives of 2 : evaluated at r = Fq. For clarity we note that 
in arriving at Eqs. (3.9) and (3.10), we expanded the radial derivatives of 2 : about Fq to 
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first-order in the eccentricity through 


*'('■) = 4 + (’■ - ’■o)4' = 4 - ‘^’■o cos(i''^)y, 


(3.11) 


where we used Eq. (3.7) in calculating the second equality; the second derivative z"{r) was 
expanded in an identical fashion. 

Notice that Eq. (3.9) specihes the angular momentum per unit mass needed at a given 
radius, characterized by fg, on a given 2D surface for elliptical or circular motion to occur. 
Equation (3.10) determines the precession parameter u for elliptical-like motion on a given 
2D surface whose surface is dehned hy z = z(r). 

The method of this subsection thus far follows closely to that of an earlier work. [6] There 
it was shown that for cylindrically symmetric surfaces with a shape prohle of the form 

z{r) oc (3.12) 


processing elliptical-like orbits are only allowed when n < 2 whereas stationary elliptical-like 
orbits can only occur for certain radii when n < 1. Interestingly when n = 1, which corre¬ 
sponds to Eq. (3.12) taking on a shape prohle of a Newtonian gravitational potential well, 
there exists no radii that will yield the stationary elliptical orbits of Newtonian gravitation. 
See Fig. 2(a) for a surface of revolution plot of a Newtonian gravitational potential well. 

In this manuscript we wish to hud the 2D cylindrically symmetric surfaces that will 
generate the Ne and GRe orbits. This corresponds to Ending the slope of the surface that 
will generate elliptical-like orbits with a constant precession parameter, u, for all radii on 
the surface. To End these surfaces we relax the condition that the radius of the orbit and the 
radial derivatives of be evaluated at r = fo and solve the difierential equation, Eq. (3.10), 
for the slope of the surface, where we treat the precession parameter u as an unspeciEed 
constant. Using this method, we obtain the general solution for the slope of the surface, 
which takes the form 


dz 

dr 


^ (i + 


- 1/2 


(3.13) 


where k is an arbitrary constant of integration and z/ is a constant precession parameter. 
Equation (3.13) can be integrated to yield the shape proEle for the 2D surface, which takes 
the form of a hypergeometric function. 

A 2D cylindrically symmetric surface with a slope given by Eq. (3.13) will generate 
elliptical-like orbits for a constant precession parameter at all radii on the surface. Notice 
that the slope of Eq. (3.13) becomes complex when u > This implies that there exists no 
cylindrically symmetric surface capable of generating elliptical-like orbits with a constant 
precession parameter greater than Using Eq. (3.8), this corresponds to a minimum 
angular separation between two apocenters of the orbital motion of ^(pmin = 27r/y^ rad ~ 
208°. 

As previously mentioned at the beginning of this section, Eq. (3.13) is left unaltered 
when the intrinsic spin of the orbiting body is included in the analysis. Thus, for a 2D 
surface whose slope is given by Eq. (3.13), both rolling and sliding objects can undergo 
elliptical-like orbital motion with a constant precession parameter at all radii on the surface 
when the right initial conditions are imposed. 



0 

-1 


1 


0 
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FIG. 1. Surface of revolution plots of z{r) vs r with the slope defined in Eq. (3.14) for k = 
—1,0,+1.[10] Each surface will generate the stationary elliptical orbits of Newtonian gravitation 
when the eccentricities of the orbits are small (Ne orbits). Eor the n = —1 surface, the slope of 
the surface diverges at r = 1. The k = 0 surface is that of an inverted cone with a slope of ^/2. 
The slope of all surfaces approach \/2 in the small r limit. It is noted that part of the n = 0, +1 
surfaces have been removed to allow for better visibility of the k = —1 surface. 

B. The 2D surfaces that generate Ne orbits 

In this subsection, we present the 2D surfaces that yield stationary elliptical orbits with 
small eccentricities (Ne orbits) for all radii on the surface. Setting v = 1 in Eq. (3.13), the 
slope of the surface takes the form 

^ = ^2(1 + (3.14) 

dr 

Notice that for negative k, r has a range 0 < r < (—I/k)^/"^, where the slope of the surface 
diverges at r = (—1 /k)^/'^. When k is set equal to zero, Eq. (3.14) reduces to an equation 
of an inverted cone with a slope of \/2. This result of an inverted cone with a slope of \/2 
giving rise to Ne orbits was previously found elsewhere. [5] See Fig. 1 for the surfaces of 
revolution with k > 0, k = 0 , and k < 0 that yield stationary elliptical-like orbits for all 
radii on its surface. 

We further note that although Eq. (3.14) yields the desired Ne orbits, these orbits do 
not, in general, obey Kepler’s third law. To hnd the corresponding Kepler-like relation for 
circular orbits on a cylindrically symmetric surface with a slope given by Eq. (3.14), we 
employ Eqs. (3.3) and (3.4), set r = r = 0 and use the fact that 0 = 27r/T for circular 

orbits. We hnd a relationship between the period and the radius of the orbit of the form 

^ (3.15) 

9 

to zeroth-order in the eccentricity. Notice that for shape prohles for the 2D surface where 
1, Eq. (3.15) takes on the approximate form 

oc r. (3.16) 
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(a)Surface of revolution plot of z{r) vs r for a 
Newtonian gravitational potential well 



(b)Surface of revolution plot of z{r) vs r with the slope 
defined in Eq. (3.14) for k = +1 


FIG. 2. (a) Surface of revolution plot of a —1/r Newtonian gravitational potential well. Although 
circular orbits on this surface obey Kepler’s third law for all radial distances, there exists no radii 
where the stationary elliptical orbits of Newtonian gravitation will occur. Hence, Kepler’s first law 
is not obeyed for orbits on this surface, (b) Surface of revolution plot with the slope defined in Eq. 
(3.14) for K = +1. This surface will generate the Ne orbits for all radial distances, obeying Kepler’s 
first law to first order in the eccentricity. Additionally, Kepler’s third law is approximately obeyed 
for these orbits when r 3> 1. Both surfaces generate particle orbits that obey Kepler’s second law 
generically. 


This corresponds to the Kepler-like relationship for an object rolling or sliding on the inside 
surface of an inverted cone. Conversely, for shape profiles for the 2D surface where 1, 

Eq. (3.15) takes on the approximate form 


n~'‘2 3 

T oc r , 


(3.17) 


which is of course Kepler’s third law for planetary motion. Notice that this i:g> 1 regime 
can only occur when k > 0 as the radius of the orbiting body resides within r < (—1 /k)^/^ 
for negative k. Hence, for an object rolling or sliding in an elliptical-like orbit on a 2D 
cylindrically symmetric surface with a slope given by Eq. (3.14), Kepler’s hrst law is obeyed 
to first-order in the eccentricity for all values of k and Kepler’s third law is obeyed only in 
the r ^ I/k^/^ regime to zeroth-order in the eccentricity. Kepler’s second law is obeyed 
generically for central-force motion to all orders in the eccentricity. See Fig. 2(b) for a 
surface of revolution plot of this surface. 


10 


























IV. PRECESSING ELLIPTICAL-LIKE ORBITS OF GENERAL RELATIVITY 
WITH SMALL ECCENTRICITIES (GRe ORBITS) 


Having found the slope of the cylindrically symmetric surface that will generate New¬ 
tonian orbits in the small eccentricity regime, we now wish to extend our study to general 
relativity (GR). Using Schwarzschild coordinates, the equations of motion for an object of 
mass m orbiting about a non-rotating, spherically symmetric object of mass M in GR are 
of the form 


e GM 
f -1- 

' O I r» 


3GMe 








= 0 


(4.1) 

(4.2) 


where, in this section, a dot indicates a derivative with respect to the proper time and i 
is the conserved angular momentum per unit mass. [7, 11, 12] It is noted that the above 
equations of motion are derivable from a relativistic Lagrangian, where here we omit the 
derivation. 

Without loss of generality, we again choose the motion of the orbiting object to reside in 
the equatorial plane as the angular momentum is conserved. Upon comparison of Eqs. (4.1) 
and (4.2) with the non-relativistic Newtonian equations of motion given by Eqs. (2.2) and 
(2.3) for a Newtonian gravitational potential, we note the presence of an additional term in 
the hrst equation of motion of Eq. (4.1). This relativistic correction term offers a small-r 
modihcation and gives rise to a precession of elliptical orbits, which will become apparent 
in what follows. 

Following a similar procedure to that outlined in Section II, Eq. (4.1) can be transformed 
into an orbital equation of motion of the form 


(fr 


dr 

dcj) 


- ^ -r + 


GM 




-r + 


3GM 


= 0 , 


(4.3) 


where we used the differential operator relation of Eq. (2.4) with the time derivative replaced 
with a derivative with respect to the proper time. We again wish to explore elliptical-like 
orbits where the eccentricities are small. We choose a perturbative solution of the form 


r{(l)) = ro(l — e cos^ucj))), 


(4.4) 


where ro and u are free parameters that are to be determined by Eq. (4.3). Inserting Eq. 
(4.4) into Eq. (4.3) and keeping terms up to hrst-order in e, we hnd that Eq. (4.4) is indeed 
a valid solution when the constants and parameters obey the relations 




V 


2 


GMro 


(i-^)“ 
V ) 


6GM 

c^ro 


(4.5) 

(4.6) 


An expansive treatment of perihelion precession for nearly circular orbits of an object resid¬ 
ing in a central potential and about a static spherically symmetric massive object in GR can 
be found elsewhere. [13] Notice that when the GR correction term GM/c^tq is set to zero, 
Eqs. (4.5) and (4.6) reduce to the results of Newtonian gravitation where = GMtq and 
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stationary elliptical orbits are recovered as z/ = 1. Also notice that for a non-zero GR cor¬ 
rection term, v only takes on values less than one. This well known result illustrates how the 
stationary elliptical orbits of Newtonian gravitation are replaced with precessing elliptical 
orbits in GR, for orbits about a non-rotating, spherically symmetric massive object. 

For a given central mass M, Eq. (4.6) implies that the deviation from stationary ellipti¬ 
cal orbits increases with decreasing radius ro, as should be expected. Hence, the innermost 
planets of our solar system experience a greater relativistic correction to the Newtonian 
orbits than the outermost planets. For the student of GR, one may notice that the rela¬ 
tion presented in Eq. (4.6) is void of a dependency on the eccentricity of the orbit. In a 
non-perturbative treatment, the precession of elliptical orbits is in fact dependent on the 
eccentricity. [11, 12] This dependency amounts to a second-order contribution for orbits with 
small eccentricities and is therefore absent in our treatment as our calculation is only valid 
to hrst-order in the eccentricity. 

Notice that when ro < 6GM/c^, Eq. (4.6) yields a complex value for v. This implies 
that elliptical-like orbits are only allowed for radii larger than this threshold radius. This 
threshold radius represents the innermost stable circular orbit (ISGO). Also notice that when 
ro < 3GM/c^, the conserved angular momentum per unit mass, i, given by Eq. (4.5) also 
becomes complex in addition to u. This implies that there are no circular orbits, stable or 
unstable, for the case of an object orbiting around a non-rotating, spherically symmetric 
central object with a radius vq < 3GM/c^. 


A. The 2D surfaces that generate GRe orbits 

In this subsection, we present the 2D surfaces that generate the GRe orbits. Adopting 
Eq. (4.6) for u and inserting this into Eq. (3.13), the slope of this surface takes the form 

I -(1 -(4.7) 

where we dehned the dimensionless relativistic quantity (5 = QGM/c^tq. The validity of this 
solution for the slope can easily be verihed by hrst inserting Eq. (3.7) into Eq. (4.7) and 
expanding to hrst-order in e and then inserting both Eqs. (3.7) and (4.7) into Eq. (3.6), 
with 1 / dehned by Eq. (4.6). 

A 2D cylindrically symmetric surface with a slope given by Eq. (4.7) will generate GRe 
orbits for all radii on the surface. Notice that the slope of this surface is dependent upon 
the central mass, M, and the characteristic radius of the orbiting body, ro, of the celestial 
system whose orbital dynamics this surface seeks to mimic via the dimensionless relativistic 
quantity 13. Further, the slope of this surface depends on [3 in both the overall factor and in 
the power of r of the 2D surface, as can be seen in Eq. (4.7). Notice that a surface dehned 
by a larger f3, which equates to a smaller value of the precession parameter u via Eq. (4.6), 
will generate an orbit with a larger precessional rate of the apocenter, as can be seen in 
Eq. (3.8). Hence, one needs a unique 2D surface to replicate each unique free particle orbit 
in GR. This result dihers dramatically from the result of the Newtonian treatment of the 
preceding section as the slope given by Eq. (3.14) is independent of the parameters of the 
theory. 

Also notice that the general relativistic relation of Eq. (4.7) reduces to its Newtonian 
counterpart given by Eq. (3.14) when the relativistic quantity (3 is set equal to zero. Al- 


12 




FIG. 3. Surface of revolution plots of z{r) vs r with slopes defined in Eq. (4.7) for k = +1 where 
the relativistic quantity /3 ranges from /? = 0 — 0.9 and increases in 0.1 intervals, displayed from 
bottom to top respectively. Each surface will generate the precessing elliptical-like orbits of general 
relativity when the eccentricities are small (GRe orbits) for the given /3. The shape of each surface, 
and hence the elliptical-like orbits each surface will generate, are dependent upon the dimensionless 
relativistic quantity [3 = 6GM/c^ro. A surface with a larger value of /3, which equates to a smaller 
value of the precession parameter u, will generate a larger precession for the orbiting body. 


though the perturbative treatment of this manuscript is only valid for orbits with small 
eccentricites, the treatment presented here is fully relativistic. The quantity f3 has a range 
0</9<l,as/9 = l corresponds to the innermost stable circular orbit (ISCO) discussed pre¬ 
viously where elliptical-like orbits are not allowed. As this relativistic quantity approaches 
one, the slope given by Eq. (4.7) diverges for all values of r. See Fig. 3 for ten stacked 
surface of revolution plots with k = +1, where the relativistic quantity (3 ranges over the 
interval (3 = 0 — 0.9, increasing in 0.1 intervals. In short, a surface dehned by a larger value 
of (3 generates an orbit with a larger precession of the apocenter. 

For the planets of our solar system, the relativistic quantity (3 has been calculated for 
each and are displayed in Table I. We note how incredibly small the relativistic parameter is 
for each of the solar system planets. These values for {3 offer a tiny modihcation to the slope 
of the surface dehned by Eq. (3.14) and replace the stationary elliptical orbits of Newtonian 
gravitation with the precessing elliptical-like orbits of GR. 


V. CONCLUSION 

In this manuscript, we showed that there exists no 2D cylindrically symmetric surface 
residing in a uniform gravitational held that is capable of reproducing the exact particle 
orbits of an object subjected to a generic 3D spherically symmetric potential. By employing 
a perturbative method to the equations of motion for an object orbiting on a 2D cylindrically 
symmetric surface, we found the general solution for the slope of the surface, given by 
Eq. (3.13), that will generate elliptical-like orbits with small eccentricities with a constant 
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Planets 

e 

ravei^) 

/? 

Mercury 

0.2056 

5.8 X 10^° 

1.5 X 10-'^ 

Venus 

0.0068 

1.1 X 10^^ 

8.3 X 10“® 

Earth 

0.0167 

1.5 X 10^^ 

6.0 X 10“® 

Mars 

0.0934 

2.3 X 10^^ 

3.9 X 10“® 

Jupiter 

0.0484 

7.8 X 10^^ 

1.2 X 10“® 

Saturn 

0.0557 

1.4 X 10^2 

6.2 X 10“^ 

Uranus 

0.0472 

2.9 X 10^2 

3.1 X 10“^ 

Neptune 

0.0086 

4.5 X 10^2 

2.0 X 10-9 


TABLE I. Table displaying the eccentricities of the orbits, the average Sun-planet distances, and 
the dimensionless relativistic quantity /3 = 6GM/c^ro for the solar system planets. [14] The small 
values of ^ for the solar system planets reveal the scale of the modifications needed to bring the 
slopes of the surfaces that generate the Ne orbits to coincide with those that generate the GRe 
orbits. 


precession parameter for all radii on the surface. We then examined this solution for the 
special cases of the surfaces that will generate Ne and GRe orbits. For a 2D surface with 
a slope given by Eq. (3.14) with k > 0, Ne orbits can be generated to hrst-order in the 
eccentricity, where Kepler’s third law will be met to lowest order in the eccentricity when 
r 3> The orbits generated on this surface differ from the elliptical-like orbits on a 

2D surface corresponding to a Newtonian gravitational potential well, where Kepler’s hrst 
law is not obeyed for any radius on the surface. For a 2D surface with a slope given by Eq. 
(4.7), GRe orbits can be generated. The slopes of these respective surfaces are functionally 
dependent upon the mass of the central object and the radius of the orbiting body of the 
gravitationally bound system whose orbital dynamics these surfaces seek to mimic. These 
surfaces reduce to their Newtonian counterpart when the relativistic term is set to zero, as 
is expected. 

A comparative study of these surfaces, whose slopes are dehned by Eqs. (3.14) and (4.7), 
offers a concrete visualization of the deviations that emerge between Newtonian gravitation 
and general relativity for an object orbiting about a spherically symmetric massive object, 
when the eccentricity is small. As embedding diagrams prove quite useful to the beginning 
student of general relativity, it can be argued that these 2D cylindrically symmetric surfaces 
that yield Ne and GRe orbits are equally useful in gaining insight into these respective 
theories of gravitation. 

An interesting undergraduate physics/engineering project could entail the manufacturing 
and experimental testing of a 2D cylindrically symmetric surface whose slope is dehned by 
Eq. (3.14). With the advent of 3D printing, the manufacturing of such a classroom-size 
surface should be feasible. Using a camera mounted above the surface and Tracker, [15] a 
video-analysis and modeling software program, the angular separation between successive 
apocenters of elliptical-like orbits with differing eccentricities could be measured and then 
compared with the theoretical prediction of 2t[. As the analysis presented in the manuscript 
is valid to hrst-order in the eccentricity, one expects a larger deviation from 27r for an orbit 
with a larger eccentricity. For details on the above experimental setup, which was used to 
analyze circular orbits on a warped spandex surface, see a previous work of the author. [4] 
For the student who seeks to generate the GRe orbits, a surface with a slope dehned by Eq. 
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(4.7) could equivalently be constructed. To generate an elliptical-like orbit that advances 
by 10° per revolution, the relativistic quantity (3 should be assigned a value of ~ 0.05 
for the slope of this surface. The inclusion of these surfaces into the classroom could offer 
the student of Newtonian gravitation and general relativity an incredibly useful tool, as 
the experimenter can see, and measure, some of the properties of these elliptical-like orbits 
hrsthand. 
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